Fragile-strong transitions and polyamorphism in glass former fluids 
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A simple model of a glass former fluid, consisting of a bidisperse mixture of penetrable spheres is 
studied. The model shows a transition from fragile to strong behavior as temperature is reduced. 
This transition is driven by the competition between the two mechanisms that contribute to diffu- 
sivity in the model: collective rearrangement of particles (responsible for the fragile behavior), and 
individual particle motion (which gives rise to the strong behavior at low temperature). We also 
observe a maximum of diffusivity as a function of pressure that can be interpreted within the same 
framework. The connection between this behavior and polyamorphism is addressed. 



I. INTRODUCTION 

The dependence of single particle diffusivity D (or al- 
ternatively viscosity 77) on temperature is the main indi- 
cator that allows to classify glass formers as "fragile" or 
"strong"l^l. For strong glass formers (such as Si02) D{T) 
is well fitted by an Arrhenius plot, indicating that there 
is a well defined value of the energy barriers which are 
relevant to diffusion. For fragile glass formers, diffusion is 
the result of collective rearrangements of the position of 
particles. As temperature is reduced, the size of the cor- 
related clusters increases- and D decreases more rapidly 
than the Arrhenius formo. There is evidence that some 
materials show a fragile behavior at high temperatures, 
whereas D{T) crosses over to a strong behavior at low 
temperatures. Water is an exampleB. The microscopic 
mechanism that originates this behavior is not well un- 
derstood. 

Some materials display the phenomenon of polyamor- 
phism, namely, the existence of different amorphous 
(solid) configurations at very low temperature. One 
well known example is, again, watercl. In this case the 
two different amorphous states transform into each other 
through a first order transition when pressure is applied 
at low temperatures. Polyamorphism has been numer- 
ically found in a model for Si02u and in some tetra- 
hedrallv coordinated materials: Si, AI2O3-Y2O3, and 
otherstl. 

One route to the appearance of polyamorphism is the 
possibility of two different local arrangement of parti- 
cles, of different density. If this is the case (and this will 
depend on the interaction potential between particles), 
external pressure can make the structure collapse from 
the most expanded configuration to the densest one, near 
some typical pressure Per- It has been recently shown in a 
simple model of spherical, purely repulsive particles with 
these characteristicaJ, that the dependence of density on 
pressure at very low temperatures has an anomalously 
rapid increase (in other words, that the isothermal com- 
pressibility has a maximum) for pressures around some 
crossover value Prr at which the two local environments 



become equally stable. We will call this a 'weak' form 
of polyamorphism. If some kind of attraction between 
particles is added, this weak polyamorphism may trans- 
form (depending on the details of the attraction) in a true 
polyamorphic first ordpr transition between two different 
amorphous structuresD. The interaction potential used in 
[7] to obtain two different local environments consists of 
a hard core at distance rg and a shoulder that extends up 
to some other distance ri > rg. Interparticle distances 
r ro are favored at high P, whereas at low P, r ri 
is preferred. 

There is an intimate relation between (weak) 
polyamorphism and fragile-strong transitions. We will 
show in this paper that a model within the class of those 
discussed in Ref. [7] , displays (weak) polyamorphism and 
has a transition from fragile to strong behavior when T is 
reduced. The model also exhibits a maximum in D as a 
function of pressure, which is well known to occur, for in- 
stance, in waterO. All these phenomena are consequences 
of the existence of two different local environments avail- 
able for the particles, and the two qualitatively different 
mechanisms that contribute to the diffusivity in these 
systems, as explained below. 



II. THE MODEL 

We studied a model similar to that of Ref [7] , with hard 
core radius vq = 0. In this way, particles are penetrable, 
and can overlap. A bidisperse system is considered, to 
avoid crystallization. Particle i is characterized by the 
value of r\ , which is taken from a bimodal distribution, 
r\ = Th, with equal probability. The 
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interaction potential between particles i and j takes the 
form 



V'^ir)=eo{l-r/R){R/ff 
V'^ (r) = for r>R 



for r < R 



(1) 



where R = r\ + r^, r is the average radius of the parti- 
cles, and So sets the energy scalecl. In all the results to be 
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presented, we use ra = 0.4375f, and rf, = 0.5625f. Tem- 
perature will be given in units of eo/ks, and pressure in 
units of eo/f^. 

We used standard Monte Carlo techniques, in the NPT 
ensemble. The system is given the dynamics of the Monte 
Carlo evolution, and diffusivities are calculated using this 
as the temporal evolution. At each time step the position 
of a single particle is modified to a new position which 
is randomly chosen within a sphere of radius 0.06 f cen- 
tered at the original position. This trial movement is 
accepted according to the Metropolis rule. The update 
is made sequentially for all particles. One Monte Carlo 
step corresponds to a sweep over all PWfifiisp. 

Previous studies of related potentiald3't2Hiil have shown 
that they display some 'anomalous' behavior, such as a 
line of maximum density and compressibility in the P-T 
diagram, and regions of anomalous melting in the P-T 
plane. These anomalies match those of real materials 
(such as water, Si02, etc.) where polyamorphism and 
fragile-strong transitions are observed or expected. We 
will concentrate here in the fluid phases of the model 
(see [10] for a discussion of the crystalline phases which 
are non-trivial) focusing on the behavior of diffusivity D. 
This is calculated from the long time behavior of the dis- 
tance traveled by the particles, < (r(t cxd) — r(O))-^ >= 
6Dt, where r{t) is the position of the particle at time t. 
The reported values always correspond to an average over 
all particles in the system. D will be given in units of 
f'^/to, where tg is the time corresponding to one Monte 
Carlo step. The number of Monte Carlo steps which are 
necessary to get good statistics greatly depends on the 
parameters P and T. As much as 5 x 10'' steps were sim- 
ulated for the lowest temperatures reported. The results 
correspond to a system of 200 particles, with periodic 
boundary conditions. 

III. RESULTS 

We start showing results indicating the existence of 
weak polyamorphism in the model. In Fig. ^ we see the 
evolution of the density p of the system and the cor- 
responding values of isothermal compressibility Kt = 
p-^dp/dP at T = and T = 0.03. Each point was 
obtained by an independent simulation at fixed pressure, 
and reducing temperature. There is a rather abrupt rais- 
ing in density and a maximum in Kt around P '--^ 1.3 
at T = 0. The maximum becomes smoother and its po- 
sition moves to slightly higher pressures at higher tem- 
peratures. The existence of this maximum is easy to un- 
derstand. At low pressures (and low temperatures) the 
system behaves as a set of hard spheres, overlapping of 
particles in the same spatial position is highly improbable 
due to energetic reasons. Upon increasing external pres- 
sure, a value Per is reached beyond which it is energeti- 
cally more favorable to accommodate particles in pairs, 
at each spatial position. The change of stability between 



'singled' and 'doubled' local configurations is ijesponsible 
for the maximum in Kt{P) around Per = 1.3E£I. This is 
consistent with the form of the radial distribution func- 
tion S{r), plotted in Fig. | for T = 0.03, at different 
values of P, where we see the rapid increase of weight 
around r ~ (indicating the appearance of 'doubled' 
particles) when crossing Per- This is the evidence for 
weak polyamorphism in the model. 
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FIG. 1. Pressure dependence of the density p and isother- 
mal compressibility Kt, at T = and T — 0.03. Lines are 
guides to the eye. 
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FIG. 2. The radial distribution function at T = 0.03, for 
three different values of P. There is a rapid increase of weight 
at distance r ~ ro = for P ~ Per- The three peaks around 
r = r originate in the bidispersity of the system. 

Next, we show results for single particle diffusivity as 
a function of pressure D{P), at different values of T in 
Fig. H We see clearly a maximum in D at P ~ Per- 
The origin of this maximum can be understood in the 
following way. For pressures much lower (larger) than 
the crossover value Per iMstly singled (doubled) config- 
uration of particles existEj. For P ^ Per, both singled 
and doubled particles coexist in the sample. Moreover, 
local rearrangements in which one particle of a doubled 
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pair moves onto some nearby singled particle are pos- 
sible. This diffusion mechanism (which is illustrated in 
Fig. H) is available only for P ^ Per, and produces an 
enhancement of diffusivity in this pressure range. 
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FIG. 3. Diffusivity as a function of pressure, at differ- 
ent temperatures. An anomalous maximum is clearly visible 
around Per. 



bination of the 'fragile' part plus this 'strong' partEj. At 
large temperature the fragile part dominates. On reduc- 
ing temperature the fragile contribution decreases much 
more rapidly, and leaves the strong contribution to be the 
main part at very low temperatures. In Fig. ^we see the 
results for D{T) in our system, at three different values 
of P. There is, in fact, a crossover between a fragile be- 
havior at high temperature and a strong behavior at low 
temperature. The fragile part of these curves can be fit- 
ted by a Vogel-Fulcher lawB-El {D - Dq exp(^/(r - To)) 
for the diffusivity of fragile glass formers, but the whole 
curve cannot be fitted to a single expression of this type, 
as we see in the figure. Also to be noticed is the fact 
that the strong component of diffusivity at low tempera- 
ture is maximized around Per- So it is close to Per than 
the fragile strong transition is more easily observed. To 
compare with a more usual case, we made an indepen- 
dent simulation for particles with a strict hard core at 
R = r\ + r^. In this case the properties of the system 
depend only on P/T (or alternatively, on density), i.e, 
the system is athermal. But we still plot D vs at 
a fixed P — 1.2 to compare with the general case. The 
results are the symbols at the left of Fig. |^. In this case 
the whole curve can be nicely fitted by a Vogel-Fulcher 
law, the fragile behavior extending down to the lowest 
values of diffusivity that we can detect. 




FIG. 4. Two-dimensional sketch of the mechanism respon- 
sible for additional diffusivity when 'singled' and 'doubled' 
particles coexist in the system. One of the particles in a pair 
(the gray one) moves from ^4 to B by successively jumping 
onto neighbor singled particles. Notice that this mechanism 
is not a collective one. 

The third set of results concerns the temperature de- 
pendence of diffusivity. In a standard glass former (i.e., 
pure hard spheres, or Lennard-Jones particles) diffusiv- 
ity occurs via rearrangements of clusters of particles. As 
temperature is reduced, the size of the typical cluster 
that can rearrange increases, producing the typical tem- 
perature dependence of diffusivity of fragile glass formers. 
This mechanism is present in our system, and is basically 
the only one active for P much larger or lower than Per- 
But for P ^ Per the mechanism discussed in the previous 
paragraph is also available. Since this is essentially a sin- 
gle particle mechanism (with a well defined activation en- 
ergy) , it account for a contribution to diffusivity which is 
basically Arrhenius-like. The total diffusivity is the com- 




FIG. 5. Diffusivity as a function of temperature for differ- 
ent pressures around Per- The high temperature part of the 
results can be fitted by a Vogel-Fulcher law, as shown for the 
P = 1.2 case. At low temperature diffusivity crosses over to a 
new regime, corresponding to a strong behavior. Left curve: 
results of a simulation for particles with a strict hard core 
at distance R, for P — 1.2. Now the whole behavior fits in 
a Volger-Fulcher law, as indicated. The circles indicates the 
points at which the snapshots of Fig. [6] were taken. 

To get a visualization of the two different mechanisms 
for diffusivity we present the snapshots of Fig. |[ They 
show the projection of the positions of the particles onto 
a reference plane, and the displacements of the particles 
during a simulation of 5 x 10^ Monte Carlo steps. Panel 
(a) was obtained at P = 1.2, T = 0.012, in our system 
of particles interacting through (|^). According to Fig. ||, 
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difFusivity has a strong behavior in this case. Panel (b) 
was obtained for the system with strict hard cores at dis- 
tance i?, at P = 1.2, r = 0.039. In this case difFusivity 
has a fragile behavior. Note that the overall values of 
D in both cases (a) and (b) are comparable (see Fig. ||). 
However, the observable features of diffusivity are clearly 
different. In the fragile case (Fig. ^(b)) a cluster of par- 
ticles (approximately delimited by the dotted line, note 
however that this is only qualitative, since the system is 
three-dimensional) has rearranged during the simulation 
time. This is consistent with the collective nature of the 
rearrangements in the fragile regime. For the strong case 
however, as Fig. ^(a) shows, we see that mostly inde- 
pendent motions of particles in uncorrelated positions of 
the sample have occur. This is the manifestation of the 
single particle character of diffusion in this case, and is 
consistent with the pictorial image of Fig. ^. 




IV. DISCUSSION AND CONCLUSIONS 

We have studied the relation between fragile-strong 
transitions and polyamorphism in a model glass former. 
We showed that if the system has a tendency to reach two 
different amorphous configurations when cooled down be- 
low and above some crossover pressure Per, then on gen- 
eral grounds we can expect maxima in the diffusivity D 
as a function of P, occurring precisely around Per- In ad- 
dition, D{T) calculated at P ^ Per is expected to have 
a fragile-strong transition when temperature is reduced. 

This behavior is in fact expected in the wide class of 
tetrahedrally coordinated materials. In all these materi- 
als, there are two typical local arrangement of particles, 
one more expanded (corresponding for instance to the lo- 
cal structure of ice Ih in the case of water) and another 
more compact (corresponding to local arrangements in 
high pressure polymorphs of ice), which correlate with 
our two typical distances between particles in the model. 
Then we expect in all these cases fragile-strong transi- 
tions and maxima of D{P) to be present. 

Real materials include in general some sort of at- 
traction between particles. With this ingredient weak 
polyamorphism may develop into a true first order phase 
transition between two different amorphous phases. This 
is in fact what happens in water, which is the best known 
example in which a truly first order polyamorphic tran- 
sition, a fragile-strong transition and maxima in D{P) 
occur. In this case the two amorphous structures are sep- 
arated by a first order line that is supposed to end in a 
critical point in the supercooled liquid region. The effects 
we have discussed (maxima in D{P) and fragile-strong 
transitions) however do not require the existence of this 
critical point. Only weak polyamorphism is needed. 
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FIG. 6. Projections of the displacements of the particles 
during a simulation of lO'^ Monte Carlo steps, at values of P 
and T as indicated. In (a) particles are taken as strict hard 
spheres of radius ri, whereas in (b) the interaction potential 
is that of (^. The difFusivity corresponds to that of a strong 
material in (a), and a fragile material in (b). 
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